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PREFACE 


The purpose of this research was to find solutions of the diffe- 
rential equation of vorticity in gradient flcw and to investigate the 
application of such solutions to the forecasting of upper air trajec- 
tcries. 

This work was conducted at the U. S. Naval Postgraduate School, 
Monterey, California, during the period December 1950 to June 195]. 

I wish to express my appreciation to Professor W. D. Duthie for 


advice and guidance throughout. 
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I. INTRODUCTION 


The possibility of making acourate mechanical forecasts of weather 
phenomena by straightforward solutions of the basic equations governing 
atmospheric motions has long intrigued meteorologists. Theoretioally, 
such a@ solution to the problem is possible but there are so many varied 
and complex forces at work in the atmosphere that all attempts so far to 
effect a complete solution to these equations have met with failure. 

One of the more important difficulties in the way of reaching satis- 
factory solutions has been the non-linearity of the equations and no 
standard methods are available to integrate nonlinear partial differential 
equations. Therefore the introduction of the so-called "perturbation" 
method was made which consists essentially in linearizing the equations by 
assuming meteorological quantities having basis values unchanging with time, 
plus small perturbation values whose second order terms oan be neglected. 

One approach to quantitative forecasting has been the linearization 
and solution of the absolute vorticity equation. This was introduced by 
Rossby [10] , who obtained a solution for the motion of sinusoidal waves 
of infinite lateral extent in a horizontal plane and found teat the velooity 


of propagation of such waves is given by the well known formula 


2 
Gz Ww, 8 4. 
4m 
where C is the velocity of the waves toward the East, U is average speed of 
the westerly current, L is wave length and # is the meridional rate of change 


of the Coriolis parameter. A considerable literature has grown out of the 
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interpretation of this formule on the classification of large scale 
atmospheric systems in terms of the so-called zonal index, which is a 
measure of the strength of the zonal current. 

Haurwitz extended Rossby's method and obtained solutions for the 
motion of waves of finite lateral extent on a horizontal plane [ 5] and 
on a sphere [6]. Craig [1] obtained solutions of the vortioity equa- 
tions in complete form without linearization for the plane and the sphere. 
These solutions differed from those of Haurwitz by the absence of a zonal 
velocity term. Neamtan [7] then made a treatment of the vorticity in 
the same manner as Craig, viz by the use of stream fwmctions and obtained 
identical solutions to those cf Haurwitz. He showed that Craig's solutions 
were inoomplete, causing absence of the zonal index term. 

Forsythe [3] developed formulas for the speed of propagation of waves 
with changing shape by use of the scalar relative vorticity as an identifiable 
property. However, he did not make any test on his formulas to determine 
their usefulness in practical forecasting. 

Rossby and co-workers [11] supplemented his initial formula based on 
the conservation of absolute vorticity with a more general technique whioh 
would be applicable to arbitrary initial streamline pattern so that the dif- 
ficulty of defining a prevailing wave length would be circumvented. 

All of the above solutions of the vorticity equations were predicated on 
the assumption of frictionless, non-divergent, and autobarotropic flow. 
Further, in order to obtain a solution for deities vorticity which could be 
useful in forecasting trajectories it was necessary for Rossby [11] to as- 


sume no horizontal shear of the wind and stationary pressure systems. 
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In this treatment, a general solution of the non-linear differential 
equation of the vertical component of absolute vorticity will be obtained 
without assuming non=-divergent flow or barotropic conditions. A solution 
is obtained for (1) a surface of constant height, and (2) a constant pres- 
sure surface. A discussion of mean-value constants for the wind shear term 
and for the effect of moving pressure systems, by means of which the vorticity 
equation could be integrated without disregarding these two terms, is included. 
& chapter is devoted to examining this solution with regard to its adaptation 
to forecasting of upperesir flow pattern by the trajectory method. In addition, 
special forms of the vorticity equation resulting from various assumptions are 
developed and discussed. 

Gradient flow is assumed throughout this investigation. The use of 
gradient flow permits a special handling of the divergence term. The wind 
of course is not always gradient even in the free atmosphere but is olocely 
approximated by the gradient wind at elevations greater than 1000 meters 
above the ground. However, under conditions of rapidly changing pressure 
gradient this close relationship between observed wind and gradient wind 
is greatly modified, and is due to the fact that the motion is not under 
balanced forces. Under these conditions there is a velocity component 
along the isallobaric gradient. Consequently the use of gradient flow in 


the vorticity equation would generally be in error. 
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II. SOLUTION OF THE VORTICITY EQUATION FOR GRADIENT FLOW 


In obtaining solutions to the differential equation of the vertical 
component of absolute vorticity as few assumptions as possible will be 
employed. Initially, only the following assumptions are made; | 

a. Friction is negligible 

be Vertical velocities are negligible. 

Other restrictions will be added later to obtain solutions of particular 
cases. The general equation of vorticity to be solved, which is due to 


Bjerknes, is: 


BS ORI 8 AY’ VA. (ed) 
play a ; 


It is noted that this equation consists of two Pale ls a term which 
yA? 


represents the horizontal divergence of the eradie-." ‘wind and a sole« 
noidal term which expresses the component of the geostrophio wind along 


the gradient of temperature. 
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be derived: 
Vv. PY = —~ VPxk wr t KV) ¢(A4FKD 7 Pxie 


which, solved for /.Y , yields 
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Now if we put this value of V.V in equation (2.2) and divide by <4 
atta ult] SA 
we obtain: DE Te a ee 


tar ae MO 
ni 
Lt. 4. AVI + Px Ie* v(n+KV)- Bre OPER) (2.3) 
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te 


‘may 


Aan a bok 
V. VPso decause © reotors, perpendicular. Hence 


we oH IT i2.. 7 yt Se 7 Px Ik) 
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Now if we make the substitution 


ba kv)v - ~~ VP«Kk 


Be 





we obtain 
1.4% iV -O+4KV V.7( Ae Ki) £ ZT. (At Kev 
2 = ae )y CT Se nde | 


Rearranging and oollecting terms 
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Integrating along the path of the particle, 
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and since V dAt= ae. 


[Aas =-(aF arty) 42. [eth (Cn+eV] 


(2.65) 


Now = 1 KVE™, where jay is the horizontal shear of the gradient 
wind, and the relation of curvature of streamlines to curvature of trajec- 
tories is KV= Keve oe « Hence equation (2.5) can be written, after 


integrating and rearranging terms, 


Se Tk +(h-KN-2% ]o exp- Sees — -K)v dln) 


Se, Toe 4 (eK V- 27, ar) 


where the subscript) indicates the initial values of the variable quantities. 
We note from equation (2.6) that on following a particle along a path 
from an initial point to some other point on the path, the ratio of their 
absolute vortioities is equal to the ratio of their temperatures times the 
inverse ratio of terms involving curvature and shear times an exponential 
term. 
Equation (2.6) was derived without assuming either zero divergence or 
a barotropic atmosphere. Furthermore, no restrictions have yet been placed 


on the shear of the wind or changing streamlines. 
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If we assume a stationary pressure systen, hae ,» as do most 


investigations of vorticity, then equation (2.6) becomes 


: t) ($3 ¥): LY 4 (7) 
= (Sa-2%) _— a K- (2 8 ) 


l. Disoussion of Orders of Magnitude of the Mean Value Constant in 
Equation (2.8) for Stationary Systems. 


If we compare the relative order of magnitude of the terms contained in 
the exponential part of equation (2.8) some conolusions can be reached re~ 
garding the value of the term J ov so that a constant oan be assigned to 
this term which will represent its average value over most mid-latitude 
paths. With suoh a oonstant the exponential terms can then be integrated. 

An examination of 700 mb. charts appears to demonstrate a radius of 
curvature of most waves of the order of 600 miles. If we take an average 
value of 20 knots for the wind velocity the order of magnitude for the 
curvature term will be 107° sec, Estimating an order of magnitude for 
the shear term in a similar manner we use for an average shear a change of 
20 knots in a distance of 600 miles. This gives an order of magnitude of 
107° seo@! for the shear term. An average value for the coriolis force 
term A would be the value at 45° latitude. At this latitude the coriolis 
term is almost exactly equal to 10™* sec™2. , Therefore, the order of magni- 
tude for the curvature and shear terms will closely approximate 107° seo 
and for the coriolis term will be 107 seom2. It is seen from this quali-~ 


tative comparison that the coriolis force is about ten times greater than 


the other two terms in influencing changes in the absolute vorticity. 
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When these values are placed in the expression ‘ 2 we obtain 
the value of -/4 when the shear term is negative and the value of VA 
when the shear is positive. 

Therefore the née; solution (2.8) can now be written as 
sa ab = 
Zils jt (x, SRY 30) (2.9) 


where Y= /q or ah ° 
tv _ eh = Mie 
The (F) term will then be (Z) or (Z) ; (24¥) differs 
[eo <7 Te 
from 2.0 by about 5%. Now, since the ratio of Z) will usually be of 
the “ay of unity the error involved in using simply ZZ T)* instead of 
(Z) will only be about 1%. 
Hence, by neglecting the exponential term represented by y* equation 


4 
(206) can be written simply as 


-(Z (L- gy) 


-s (L-2) (2.10) 


2-e Discussion of Constants in the Equation for Moving Pressure Systems. 
All treatments of the vorticity equation appearing in the literature 
so far have ignored the effects of difference between the curvatures of 
streamlines and trajectories a » and have assumed a steady state con- 
dition. This implies stationary pressure systems which of course is not 
the prevailing case in nature. The ‘difficulty of handling gy mathe- 


matically has been responsible for neglecting this term. 
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Perhaps some conclusions can be drawn, however, from an examination 
of the range of variation of v(K-K) so that a mean value constant 
can be assigned to € with the result that equation (27) can be inte- 
grated. 

From the relationship of trajectory,ourvature and streamline curvature 
according to Petterssen [8, p. 225] , Kt Ler S tes) ) 


we can solve for V(iK-K) so that 


V(K- Ks) = Ks Clos $ 


Consider now a wave shaped system which is more or less symmetrical 


about a latitude circle and moving eastward with speed C. Along the trough 





which is to the right of the path of the system Cos ¥7° and is unity at 
bottom of troughe Along the ridge where kK is anticyclonic to left of 
path Cos %<o. If we choose average values of fo for Cos /, the sign 
depending on whether along a trough or ridge, and take average values of 
1000 km and 20 kts for &% and C respectively then Y(K-K) ~ 10” eindiye 
of magnitude. 

From previous considerations of orders of magnitude for the shear and 
absolute vorticity we arrived at values of 107° seo7! and 10™ seont 
respectively for two terms. Therefore the final order of magnitude for 


will be two possible valves: 
— +¥ me 
gif (2 —/0 
E ru 70-4 
If the shear is negative, then €7~-0.2  , and if shear is positive, 


EO, 
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Finally, substituting these two possible values in equation (2.7) 





we obtain 
2.2 f 
kus 4 2S. ), } CNp-va0: B 
a “ 1.2 Ja) 
and Se (Fy 3. PaO 
‘ e6)< yon 


Therefore, under the conditions stated above the existence of positive 
or cyclonic wind shear cancels out the effect of the a term and with 
anticyclonic or negative shear the correction term € is significant. The 
opposite effects would occur in the case of retrograde waves or in those 


parts of closed pressure systems where the wind blows in a direction 


opposite to the direction of motion of the system. 


(10 ) 
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III. SPECIAL CASES OF THE GENERAL SOLUTION OF THE 
VORTICITY EQUATION FOR STATIONARY SYSTEMS 


Certain simplifying assumptions can now be made with regard to the 


various terms of the general solution (2.8). The special cases which will 


be discussed here are: 
1. Geostrophic Flow 
2. Wind Flow Parallel to Isotherms 
Se No Horizontal Wind Shear 
4. Combinations of the Above Three Cases. 
le Case of Geostrophic Flow. 
straight isobars or contours, 


Geostrophio flow prevails where there exist. 


hence the ourvature term becomes zero and equation (2.8) becomes: 


(361) 


However, as demonstrated in Chapter II, the exponential term can be 


neglected with only a small error involved so that for the geostrophic flow 


case equation (3.1) can be written as 


zie z T* Ae 
YH, TR 5 (3.2) 





Equation (3.2) can be written as a quadratic and solved by the quadratio 


formula so that 


(343) 


(12) 
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Equation (3.3) can now be expanded binomially: 
Sf ps2 | 
Se = L641 $244 (8k) Bley (4) (Z )e*4 
-1/ 6 
+ fe (66) ©) + -—-J 


All terms beyond the first in this binomial expansion are of the order 
of 107° or smaller and which, because of a much smaller order of magnitude 
than for the shear term, oan be neglected so that for the case of geostrophic 
flow the absolute vorticity can be expressed simply by a Barty, o> Ov . 

Inasmuch as the absolute vorticity cannot be zero except for the very 
improbable case of geostrophic flow with no shear and at the equator, the 
choice of sign on the second term in (3.3) must be positive. 

For practical application on the weather chart in determing how such 
flow would be affeoted by various types of movement, the relative vorticity 


is examined. For geostrophic flow, (3.2) can be written 


“ i ) 
sophie 


A particle moving northward would therefore undergo a decrease in 
relative vorticity or be turned anticycolonically even though there may be 
a slight positive contribution from the temperature term. Southward flow 
would produce the opposite effeot. 

A particle moving along a parallel of latitude would be influenced only 
by the change in temperature. If moving toward higher temperature, there 
will be an increase in relative vorticity and particle will ourve to the 


left. 


(12 ) 
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2e Case of Wind Flow Parallel to Isotherms. 

The wind is blowing parallel to the isotherms in this oase, which 
implies a barotropio atmosphere. Such conditions prevail when systems 
are thermally symmetrical. In this case equation (2.8) becomes 


te - (RV eae 


a. (KY +42) (3.4) 


and the exponential term becomes unitye Now, on following a particle from 
initial point or path to some other point the pcluinice vorticity changes 
as the inverse ratio of their respective (Kv #0) terms e 

If we now expand equation (3.4) a quadratic is obtained which can be 
written as y* Joe la » where c=(Kevtr 4a) (Lv4>)- , the 
initial point values. Solving this equation by means of the quadratic 
formula gives 


Re GANA en da. 


(3.5) 


Equation (3.5) can now be expanded in a binomial expansion giving in a 


manner similar to Case ls 


—) 
= (548) 


The choice of the positive sign for the second term is discussed under 
Case l- 


Equation (3.4) can be now written as relative vorticity 


Ji (Keven), KV40 422). 7 
Ken). 


(13) 













pee A ALpeet weet ore eee. 


doe «9240 @08e ht Serer te) ee Lee peel ot eee ee 


iy, Firs, oma = se Lo 

Sie i eo he ee 
‘Ysug™- ?. 
ame bev = 


et jeri 1 li) © ee a Re 2 . ieee ae i’ 
cers ph feds WOiwets “ek! Tey et Of Arey om pace Eten 
dlink Ne BR) oer Wah ape ee eee 
en i) image Goer eee 
x Raat omer. 5a “2 be eso 
ON St Ree fm ee tweeted ormney trie Lena 
—\ 


G =e: Bw oo 



















ee er 
“eee eka 


4 
=< "thle iain di i ee te eee 
—— eee 


See peees * 











ad 








rath 





- 


3e Case of No Horizontal Wind Shear. 
If the horizontal shear of the wind is either assumed to be zero or 
is neglected then the general equation (2.8) simplifies to 


Se. T (Keven), 
Tt. Te” (KiV40) (3.7) 


This oan be written as 


thant. a8 ff 





= 
—_ 


a ! 
oy ns 2 ——= 8 jec 
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We see then that for the oase of no shear the ratio of the absolute 
vorticity at two points is equal to the ratio of their respective absolute 
temperatures. It is also to be noted that this relationship whioh involves 
no assumptions as to divergence of the wind or barotropy differs from Rossby's 
formula as applied to the forecasting of particle trajectories only by the 
presence of the temperature terme 

If we express the vorticity as ae and choose an inflection latitude, 
as did Rossby and oo-workers [11] , at which the ourvature is zero, then 
equation (3.8) oan be expressed as 

77 


KsV= Ne =D (3.9) 


which corresponds to Rossby's formula 


a TI : 


Equation (3.9) could be applied in a manner entirely analogous to the 
technique used by Rossby in forecasting air partiole trajectories. The 


coriolis force at any point from the inflection latitude can be expressed 


(14) 
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as a function of distance from such latitude so that d= By + De where & 
is the rate of change of coriolis force with latitude and is considered 


constant over the latitude range concerned. Then (3.9) can be written as 


eve E de - (A, +2.) ‘ (3210) 


If we express the curvature in terms of a second order differential 
equation of its Y coordinate or distance from the inflection latitude we 


obtain after substituting this in (3.10) 


hee Bo 
Te ne (Z-1) -4, 


This differential equation could perhaps be solved by elliptic integrals 
if a proper constant is chosen for the temperature term so that a value of y 
coordinate would be obtained as a function of the intersection angle of 
particle path with inflection latitude. 

However, a qualitative interpretation of equation (3.10) can be made 
showing effects on the relative vorticity with varying paths without actually 
solving the differential equation. If the particle is moving northward init- 
jally, the effect of coriolis force is to decrease the relative vorticity 
and inoreasing temperature tends also to increase the vorticity. However, 
now assuming decrease in temperature northward the temperature term in 
general is of much smaller magnitude el , hence the particle curves 
anticyolonically but to a slightly greater extent than the path resulting 
from coriolis considerations only. If the temperature increases northward 
along the path the effect is to decrease the effect of inoreasing coriolis 


force so that relative vorticity is decreased to a lesser extent than when 


(15) 
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the coriolis effect is considered alone. Therefore the path would curve 
less anticyclonioally than otherwise. 

For a southward moving current the converse of all the above state- 
ments would hold. 
4. Combinations of the Three Preceding Cases. 

a. Geostrophic Flow with no Wind Shear. 


Under such conditions of wind flow equation (2.8) becomes simply 


ta T* Az 
Sy FP 
which can be simplified to 
te a s a 
a ee 


or the change of met Srey ck is now due simply to change of 
coriolis force)is directly proportional to the ohange of respective 
absolute temperatures. It is also interesting to note that the flow 
in this case implies a gradient of temperature from South to North, 
which of course, is not true in nature except under local conditions. 
be Geostrophic Flow with a Barotropic Atmosphere. 


In this case the general equation (2.8) can be written 


ae a 
. Re odes os 


V 
(n+Bt). 2 
Ce A Barotropio Atmosphere with no Wind Shear. 


Under these conditions equation (2.8) becomes 
> Ks V+ >). 


Ja, (Kevan) 


This can be written as 


(ever =Kevire . 
(16) 
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This result is identical in form to that used by Rossby and co- 
workers [11] in foreoasting constant vortioity trajectories with the 
significant difference that no restrictions have been placed on divergence 
or convergence in this development. Therefore, the same technique as used 
by Rossby in forecasting particle trajectories could be applied to the 


gradient wind using identical relationships as: 


Kv >-F7 


where = rate of change of coriolis force and (y) equals distance from 
an inflection latitude at which relative vorticity and curvature are zero. 
Solutions of the second order differential equation for the maximum y co- 
ordinate as a function of the intersection angle of the path with the in- 
flection latitude would be identical to those of Rossby [11] or Fultz [4]. 
The verification of forecasts of such computed trajectories should be 
more accurate than that of Rossby because of the added refinement of no 
restriction on divergence and convergence. 


d. Combination of all Three Special Cases -- Geostrophio Flow, No 
Wind Shear, and Barotropy. 


These conditions would give a form for the vorticity equation thus 
4. fe 


—— oe 
4, 7 
which gives the rather trivial result 
However, absolute vorticity is again shown to be constant although only 


a function of latitude. Furthermore, zonal flow is implied by this 


equality of coriolis terms. 
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IV. SOLUTION OF THE VORTICITY EQUATION FOR A CONSTANT 
PRESSURE SURFACE 
The equation of absolute vorticity along a surface of constant 
pressure is somewhat more simplified than along a oonstant height sur- 
face because the solenoid term is absent. This is not surprising since 
an isobaric surface obviously oannot be intersected by solenoids. 
Thus, equation (2.2) when applied to an isobaric surface involves 


only a divergence term: 


dk _ Kay 
At (401) 
We can now express this divergence in terms of ooriolis force, cur- 
vature and velocity in a manner entirely analogous to that for a constant 
height surface. 
Starting with the expression for the gradient wind: 
Ve Vys- (viv 


which we can express for an isobaric surface as 


V=-94 V¥xR- kv) ¥ 


(402) 
a 


Now multiplying by 7 and rearranging: 
~9V,8xk=2AV+ KVY 


solving for VV gives 
“/ 
V 9 Wee xR (M4KY) 
Taking the divergence of both sides gives 


VV =~ 9 Vk vo(2 #kv) - Ot BY) 9 Vp Bx 


(28 ) 
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which reduces to 


Ve = -9 FR VY (atkv) 


(4.3) 


sinoe Vy Vy Xx h=a0. 


It is seen that equation (4-3) gives the divergence of the gradient 
wind as the scalar product of a vector parallel to the isobypses and the 
-! 
gradient of the quantity (A+ ev) ; 


Upon substituting this +7 phen for the eg te in (4.1) and 


a ai n 
oy SG Lam Labay ae! 


ri 


rearranging we get 


+ mae + 9% Bank. Yaa lrentil 


Since ~ FER = 4V (At Kv) we can write, after performing 


indicated differential operations: 


AS. 4y. v,(atKv) 
© an Cate) 


which can be expressed as the integral 


[4Q.5) » + [LAr 


(4.4) 
Now eince AP _V and Va: Avs Aa 
at 
we oan write (4.4) as 
fae Se) = + [ch bu(ntv) 
which is readily integrable to 
AK. be Cote. 
mT (at Ko 
and =o 4 (n+ Gv) (4.5) 
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Equation (4.5) is identical with (3.4) which was obtained for a 
constant level surface by assuming a barotropic atmosphere. However, in 
developing equation (4.5) no such assumption was necessary. For equation 
(4.5) to be applicable, however, the flow must be assumed to take place 
along an isobaric surface. Therefore its application would be highly 
limited because of the improbability of isobaric flow. 

Now equation (4.5) can be rearranged as a quadratic and expanded 
binomially in a manner similar to equation (3.4) and which gives identical 


OA 


results: 


As a suggested extension to this work solutions of the vorticity 
equation could be carried out in an entirely analogous manner for isentropic 
surfaces and would have wider application than the solution for a constant 
pressure surface because of the greater prevalence of flow along isentropic 


surfaces. 


(20) 
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